Applied Energy 126 (2014) 78-89 



gataiHBBr 

ELSEVIER 


Contents lists available at ScienceDirect 

Applied Energy 

journal homepage: www.elsevier.com/locate/apenergy 



A novel thermomechanical energy conversion cycle 

Ian M. McKinley, Felix Y. Lee, Laurent Pilon* 



CrossMark 


Mechanical and Aerospace Engineering Department, Henry Samueli School of Engineering and Applied Science, University of California, Los Angeles, Los Angeles, CA 90095, USA 


HIGHLIGHTS 


• Demonstration of a novel cycle converting thermal and mechanical energy directly into electrical energy. 

• The new cycle is adaptable to changing thermal and mechanical conditions. 

• The new cycle can generate electrical power at temperatures below those of other pyroelectric power cycles. 

• The new cycle can generate larger electrical power than traditional mechanical cycles using piezoelectric materials. 
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This paper presents a new power cycle for direct conversion of thermomechanical energy into electrical 
energy performed on pyroelectric materials. It consists sequentially of (i) an isothermal electric poling 
process performed under zero stress followed by (ii) a combined uniaxial compressive stress and heating 
process, (iii) an isothermal electric de-poling process under uniaxial stress, and finally (iv) the removal of 
compressive stress during a cooling process. The new cycle was demonstrated experimentally on 
[001]-poled PMN-28PT single crystals. The maximum power and energy densities obtained were 
41 W/L and 41 J/L/cycle respectively for cold and hot source temperatures of 22 and 130 °C, electric field 
between 0.2 and 0.95 MV/m, and with uniaxial load of 35.56 MPa at frequency of 1 Hz. The performance 
and constraints on the operating conditions of the new cycle were compared with those of the Olsen 
cycle. The new cycle was able to generate power at temperatures below those of the Olsen cycle. In 
addition, the new power cycle can adapt to changing thermal and mechanical conditions. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Harvesting thermal and/or mechanical energies that are typi¬ 
cally wasted could contribute to more sustainable and efficient 
energy use. Waste mechanical energy is an unavoidable by-prod¬ 
uct of objects in motion and exists in the form of vibrations, shocks, 
or strains [1], Sources of waste mechanical energy include fluid 
flow, household appliances, industrial equipment, motor vehicles, 
and structures such as buildings and bridges [1 ]. In addition, waste 
heat is the inevitable by-product of power, refrigeration, and heat 
pump cycles, according to the second law of thermodynamics [2], 
In fact, many sources of waste heat, such as electricity generation 
and transportation systems, also waste mechanical energy. 

The most widely used method to harvest mechanical energy is 
piezoelectric energy conversion [3]. It makes use of the piezoelectric 
effect to convert time-dependent mechanical deformations into 
electricity [1], Other methods for direct mechanical to electrical 
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energy conversion include electromagnet, electrostatic, and electro¬ 
active polymer generators [1]. Similarly, various methods are avail¬ 
able to harvest waste heat. For example, Stirling engines [4] and 
organic Rankine cycles [5] have been used to convert low grade 
thermal energy into mechanical energy. Thermoelectric devices 
convert a steady-state temperature difference at the junction of 
two dissimilar metals or semiconductors into electrical energy [6], 
By contrast, the Olsen cycle [7] performed on pyroelectric materials 
utilizes time-dependent temperature oscillations to convert ther¬ 
mal energy directly into electricity. Note that none of these energy 
conversion methods are capable of converting both thermal and 
mechanical energies directly into electricity. 

Pyroelectric materials possess a temperature-dependent spon¬ 
taneous polarization, defined as the average electric dipole 
moment per unit volume, in absence of an applied electric field 
[8], A subclass of pyroelectric materials, known as ferroelectrics, 
has the ability to switch the direction and magnitude of the spon¬ 
taneous polarization by reversing the applied electric field above 
the coercive electric field [[9]. Pyroelectric materials are also piezo¬ 
electric, i.e., the electric charge at the material surface changes 
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when the material is mechanically deformed. Applying a compres¬ 
sive stress in the poling direction decreases the electric displace¬ 
ment for a given temperature and electric field. The reduction in 
surface charge results in current in the external load. 

This study aims to demonstrate the feasibility of a novel power 
cycle performed on pyroelectric materials to convert both thermal 
and mechanical energies directly into electrical energy. This new 
cycle was demonstrated on commercially available [001]-poled 
lead magnesium niobate-lead titanate 0.72PbMg 1 /3Nb 2 /30 3 _ 
-0.28PbTi0 3 (PMN-28PT) single crystals. 

2. Background 

2.3. Dielectric hysteresis loops 

Fig. la and b show the isothermal bipolar hysteresis curves in 
the electric displacement D vs electric field E diagram (D-E loops) 
exhibited by a typical pyroelectric material at two different tem¬ 
peratures T cold and T hol under compressive stress a equal to 0 and 
cr H . The D-E loops under any compressive stress traveled in a 
counter-clockwise direction. The electric displacement D of a 
pyroelectric material at temperature T under electric field E and 
compressive stress a can be expressed as [9,10] 

D(E,T, a) = e 0 E r (T, cr)E + P s (T,cr) (1) 


where e 0 is the vacuum permittivity (=8.854 x 1CT 12 F/m) and 
£ r (T, a) is the large-field relative permittivity of the material at tem¬ 
perature T and under stress a. The saturation polarization P S (T, a) 
corresponds to the electric displacement in the linear fit of D vs E 
at large field extrapolated to zero electric field 11] and the slope 
of this linear fit is the product £ 0 £ r (T, a). Other important properties 
include (i) the remnant polarization P r (T, a) corresponding to the 
polarization under zero applied electric field, (ii) the coercive field 
£ C (T, a) corresponding to the electric field required to reach zero 
electric displacement, and (iii) the Curie temperature T Cu ne defined 
as the temperature at which a ferroelectric material undergoes a 
phase transition from ferroelectric to paraelectric. This phase 
transition temperature is typically defined as the temperature corre¬ 
sponding to the maximum of the real part of the complex dielectric 
constant for given frequency and applied electric field [12], 

2.2. Single crystal PMN-PT 

Single crystal PMN-xPT has been widely used in mechanical 
sensors and actuators and their piezoelectric and dielectric proper¬ 
ties have been studied extensively 11,13-25], PMN-xPT possesses 
large piezoelectric constants near the morphotropic phase bound¬ 
ary (MPB) separating rhombohedral and tetragonal phases 20], 
This phase boundary in PMN-xPT corresponds to x ranging 
between 27.5 and 33 mol% [19], 
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Fig. 1 . Two-dimensional projections of (a) the Olsen cycle and of (b) the new power cycle in the D-E plane and electric displacement vs electric field loops at T co id for uniaxial 
stress a = 0 and T hot for a = 0 and o H . The electrical energy generated per cycle is represented by the grey areas enclosed by states 1-2-3-4. (c) The thermal, electrical, and 
stress states of the ferroelectric sample at each state of the new power cycle. 
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Properties of PMN-xPT with x between 30 and 33 mol% as func¬ 
tions of electric field and mechanical stress have been recently 
investigated! 13,1 7], For example, Feng et al. [13] reported the 
effect of uniaxial compressive stress <r, applied in the poling 
direction, on (i) the bipolar D-E loops at room temperature, (ii) 
the coercive electric field E c , and (iii) the remnant polarization P r 
for [001 [-oriented PMN-30PT single crystals [13], The remnant 
polarization P r was approximately 0.25 C/m 2 under zero loading 
but vanished for bias compressive stress a around 50 MPa. 
Similarly, Wan et al. [17] reported remnant polarizations of 
0.25 and 0.30 C/m 2 at room temperature under zero load for 
PMN-32PT and PMN-33PT, respectively. Under compressive stress 
larger than 30 MPa, the remnant polarization was negligibly small 
(~0.08 C/m 2 ). Overall, the coercive field E c and remnant polariza¬ 
tion P r decreased with increasing compressive stress [17,13], 

Recent studies demonstrated that not only external electric 
field but also mechanical loading contribute to polarization switch¬ 
ing and phase transitions in PMN-xPT single crystals [13,14,17,18], 
Zhou et al. [23] compared the performance of PMN-xPT with x 
equal to 28%, 30%, and 32% and determined that PMN-28PT was 
the preferred composition for transducer, sensor, and actuator 
applications due to its larger piezoelectric properties in a broad 
temperature range. Based on these observations, single crystal 
PMN-28PT samples were used in the present study. The electric 
field vs temperature ( E-T ) phase diagram for [001] PMN-28PT 
tinder zero stress indicates that the material can assume four dif¬ 
ferent crystalline phases in the temperature range from 27 to 
177 °C for electric field between 0 and 1.5 MV/m including the 
monoclinic phases M A and M c , the tetragonal (T), and the cubic 
phase (C) [22]. The associated phases boundaries are both temper¬ 
ature and electric field dependent. At room temperature under 
zero electric field, the material is in the rhombohedral phase 
[26], The Curie temperature at zero field was reported to be around 
125 °C [22] above which a pseudocubic state existed featuring local 
lattice distortions [27]. An applied electric field can cause a phase 
transition from the pseudocubic phase to the ferroelectric state 
resulting in double hysteresis loops [28], 

2.3. Thermomechanical cycles 

Traditional thermodynamic power cycles, such as the Stirling, 
Rankine, and Ericsson cycles, are performed on a working fluid 
and convert thermal energy into mechanical work. Their efficiency 
is defined as the ratio of the net mechanical work produced by the 
cycle to the thermal energy input. These cycles are often repre¬ 
sented in the pressure vs specific volume ( p-v ) diagrams. In some 
of them (e.g., the Rankine cycle), the working fluid undergoes 
liquid to vapor phase transition. Additionally, thermoelectric 
devices convert thermal energy to electrical energy by taking 
advantage of the Seebeck effect [6[. Their conversion efficiency is 
defined as the ratio of the electrical energy generated to the ther¬ 
mal energy provided at the hot junction [6], On the other hand, in 
piezoelectric energy conversion, mechanical energy is converted 
into electrical energy. Then, the conversion efficiency is defined 
as the ratio of the generated electrical energy to the applied 
mechanical energy [29], 

Alternatively, the Olsen cycle can be performed on pyroelectric 
materials and consists of two isothermal and two isoelectric field 
processes in the displacement versus electric field diagram [7[. It 
is analogous to the Ericsson cycle defined in the p-v diagram of a 
working fluid [7[. The first process of the Olsen cycle corresponds 
to an increase in electric field from E L to E H at constant temperature 
Tcoid- The second process consists of heating the material from T cM 
to T/, or under constant electric field E H . The third process corre¬ 
sponds to a decrease in the electric field from E H to E L at constant 
temperature T hot . Finally, the cycle is closed by cooling the material 


from Thor to T coU under constant electric field E L . Note that these 
processes are performed under zero mechanical stress. 

Fig. la shows the Olsen cycle in the two-dimensional D-E 
diagram for <7 = 0 MPa. The area enclosed by the four processes pre¬ 
viously described corresponds to the electrical energy generated 
per unit volume of material per cycle, denoted by N D . It is expressed 
in J/L/cycle (1 J/L/cycle = 1 kj/m 3 /cycle) and defined as [7] 

N d = j EdD (2) 

The power density P D (in W/L) is the amount of electrical energy 
generated by the pyroelectric material per unit volume of material 
per unit time. It is expressed as P D = N D f where / is the cycle 
frequency. 

In order to maximize the energy density N D generated by the 
Olsen cycle, the electric field span (E H - E L ) should be as large as 
possible without causing de-poling and/or electric breakdown 
[30]. Similarly, the electric displacement span should be the largest 
possible. To do so, T hot should be greater than the Curie temperature 
E curie- Indeed, when the material is paraelectric, the spontaneous 
polarization vanishes. However, near T Cur ie an applied electric field 
can re-pole the material. As the temperature increases above T Cur ie . 
the paraelectric state becomes more and more stable and the elec¬ 
tric field required to induce the ferroelectric state increases [22], 
Thus, the cycle should be performed with T hot sufficiently above 
Tcurie so that even a high electric field cannot re-pole the material. 

Numerous studies have investigated the Olsen cycle performed 
on pyroelectric polymers [31,39-41], single crystals [32,33,42] and 
ceramics [7,30,34,35,36-38], In particular, Kandilian et al. [32] 
studied the pyroelectric energy generation of 140 pm thick 
[001 [-oriented PMN-32PT single crystals using the Olsen cycle. 
The maximum generated energy density was 100 J/L/cycle for elec¬ 
tric field cycled between 0.2 and 0.9 MV/m and temperatures vary¬ 
ing between 80 and 170 °C. In addition, T Cu ne under zero field was 
determined to be around 150 °C [32]. Moreover, the PMN-32PT 
samples experienced dielectric breakdown for electric fields larger 
than ~0.9 MV/m, and samples suffered from cracking due to ther¬ 
mal stress for temperature differences in excess of 90 °C [32], 

The use of heat conduction to heat and cool a pyroelectric mate¬ 
rial in the Olsen cycle was explored by performing a “stamping 
procedure” [41]. It consisted of alternately pressing a pyroelectric 
sample in thermal contact with cold and hot aluminum blocks 
under specified electric fields. Lee et al. [41] used co-polymer 
60/40 P (VDF-TrFE) films as the pyroelectric material and achieved 
a maximum energy density of 155 J/L/cycle at 0.066 Hz for temper¬ 
atures between 25 and 110 °C electric fields between 20 and 
35 MV/m. Note that single crystals and ceramics offer advantages 
over polymers in that (i) they possess significantly lower leakage 
current thanks to their higher electrical resistivity and (ii) they 
do not require electrical poling prior to performing the Olsen cycle. 
However, polymer films can sustain significantly larger electric 
field. 

The Olsen cycle has been demonstrated to generate the largest 
energy densities of any power cycle performed on ferroelectric 
materials [43], However, the Olsen cycle performed on a given 
material requires the hot source temperature to be larger than 
Tcurie in order to generate energy. Moreover, the power density of 
the Olsen cycle is limited by the low cycle frequency due to slow 
heat transfer and relaxation processes. The energy and power den¬ 
sities generated by a pyroelectric material undergoing the Olsen 
cycle could be increased by increasing the electric displacement 
span between temperatures T coM and T hol . This can be achieved 
by subjecting the sample to compressive stress. This is particularly 
true at temperatures around T Cu rie- 

Here, we present a new cycle using uniaxial compressive stress 
in addition to thermal and electric field cycling. It circumvents the 
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above mentioned difficulties to increase both energy and power 
densities (i) by combining piezoelectric and pyroelectric energy 
conversion and (ii) by increasing the cycle frequency by quickly 
forcing the material into a specific phase using mechanical stress 
instead of heating and cooling which are inherently slow. Single 
crystal PMN-28PT was chosen for its advantageous piezoelectric 
properties. 

3. New thermomechanical power cycle 

Fig. lb illustrates the new power cycle projected onto the D-E 
plane overlaid with the corresponding isothermal bipolar hystere¬ 
sis curves ( D-E loops) at cold temperature T coW under zero stress 
and at hot temperature T hot under compressive stress a H . Fig. lc 
schematically illustrates the practical implementation of the new 
power cycle and the stress, temperature, and electric field imposed 
in each state. Process 1-2 consists of an isothermal increase in elec¬ 
tric field from E L to E H at T roW in the absence of compressive stress. 
Process 2-3 corresponds to simultaneously compressing the sam¬ 
ple at a H and heating it up to T hot . Process 3-4 consists of an iso¬ 
thermal decrease in electric field from E H to E L at T hot under 
compressive bias stress a H . Finally, process 4-1 closes the cycle 
by simultaneously cooling the sample to T coW and removing the 
loading under constant electric field E,. The area enclosed by the 
four processes in the D-E diagram, shown in Fig. lb, corresponds 
to the generated energy density N D defined by Eq. (2). The overall 
cycle frequency (in Hz) is defined as /= (T 12 ± T 23 ± t 34 ± t 4 i)~’ 
with Ty corresponding to the duration of process i-j. This new cycle 
can be implemented using a procedure similar to the “stamping 
procedure” [41 j previously discussed but under significantly larger 
compressive stress. This was implemented experimentally to 
directly convert both thermal and mechanical energies into electri¬ 
cal energy using single crystal PMN-28PT. 

3.2. Material efficiency 

The material efficiency of a power cycle is typically defined as 
the ratio of the electrical energy generated by the material during 
the cycle to the thermal and/or mechanical energy consumed by 
the material to be converted by the cycle. In the present cycle, both 
thermal and mechanical energies are converted into electrical 
energy. Then, the material efficiency of the new cycle can be 
expressed as 


Note that this definition of efficiency is consistent with that used for 
piezoelectric materials given by 17 = N D /W in [29] and for the Olsen 
cycle given by i] = JV D /Q,„ [44]. The thermal energy provided to the 
material during the cycle may be expressed, per unit volume, as 

Qin = f pc P (T)dT (4) 

where p and c p (T) are the density and specific heat of the ferroelec¬ 
tric material in kg/m 3 and J/kg K, respectively. The specific heat of 
ferroelectric materials c p (T) is temperature-dependent and can be 
measured by differential scanning calorimetry [45], The net 
mechanical work W in provided during the cycle per unit volume 
of material can be defined as 

W, n = f <x(e)de (5) 

where e represents the strain in the longitudinal direction parallel to 
the polarization. For ferroelectric materials undergoing phase transi- 
tion(s), the relationship between a and e is typically non-linear 


[13,46], Thus, W in cannot be expressed in terms of the Young’s mod¬ 
ulus. Instead, it should be estimated from stress-strain curves [47], 
Note that the material efficiency defined above accounts for the con¬ 
version of thermomechanical energy into electricity by the material 
itself. It does not represent the efficiency of a potential device imple¬ 
menting the cycle and subject to heat losses, friction, and other irre¬ 
versible processes. In other words, t] represents the upper limit of a 
device’s efficiency. The latter will also depend on the size of the 
system. 

4. Materials and methods 

4.1. Samples 

Single crystal samples of PMN-28PT were purchased from 
Sinoceramics, LLC. The dimensions of the samples were 
5x5x3 mm 3 . The average weight of the samples was 
588.5 ± 0.8 mg, corresponding to a density of p = 7847 ±11 kg/m 3 . 
These samples were poled in the [001 [-direction. The two 
5x5 mm 2 faces of each sample were coated with Cr/Au electrodes. 

4.2. Experimental setup 

The experimental setup consisted of an electrical and a 
thermomechanical subsystem. The electrical subsystem was a 
Sawyer-Tower circuit identical to that used in our previous studies 
[30-33,42,48], Fig. 2a and b show a schematic and a photograph of 
the thermomechanical subsystem used to perform the novel power 
cycle, respectively. This subsystem consisted of a spring return air 
cylinder (McMaster-Carr 6498K252) vertically actuated using 
compressed air at a maximum pressure of 469 kPa. The extension 
and contraction of the cylinder rod were controlled by a 24 V DC 
solenoid valve. The heat source consisted of a 100-W cartridge hea¬ 
ter imbedded in a 1.27 cm thick aluminum plate. A type-K thermo¬ 
couple was placed at the center of the heating block. The block’s 
temperature was kept constant using a proportional integral 
derivative (PID) temperature controller (Omega CN-7823). The 
PMN-28PT sample was sandwiched between two copper tapes to 
provide electrical contact between the sample’s electrodes and 
the wires. This assembly was placed on top of a 5 mm thick steel 
die. An aluminum heat sink (Cool Innovations 3-151514 M) was 
placed in thermal contact with the steel die by epoxy adhesive 
OMEGABOND* 200 to passively cool the pyroelectric sample to 
T cold during process 4-1. A 140 pm thick Kapton film was used to 
electrically isolate the sample’s electrodes from the metallic heat 
source and sink. Note that the sample temperature could not be 
measured during electric field cycling due to electrical conduction 
between the sample and the thermocouple. 

This setup was modified slightly when collecting isothermal 
D-E loops under different compressive stresses. Then, the heat sink 
and steel die were replaced with an aluminum heating plate and 
wood block identical to the heat source placed above the sample. 
This was done to ensure the sample was uniformly heated to T h „, 
from both sides and to minimize the temperature gradient in the 
sample. In this case, a type-K thermocouple placed directly on 
the sample was used to measure the sample temperature. It was 
removed before electric field cycling began. 

4.3. Experimental procedure 

Isothermal D-E loops under different compressive stresses were 
collected and the new power cycle was performed on the PMN- 
28PT samples using the above-described experimental setup. For 
comparison purposes, the Olsen cycle was also performed using 



82 


l.M. McKinley et al./Applied Energy 126 (2014) 78-89 




Solenoid 

valve 

control 


Fig. 2. (a) Schematic and (b) photograph of the thermomechanical subsystem used 
to create the periodic temperature oscillations and to apply compressive stress 
during the new thermoelectromechanical power cycle depicted in Fig. lb. Dimen¬ 
sions in (a) are not to scale. 


aluminum sample. The specific heat c p of aluminum was found to 
be 898 J/kg K at 27 °C falling within 0.6% of the value reported in 
the literature [50]. 

4.3.3. New power cycle 

The new power cycle was performed so that the duration of 
each process was t 2 3 = t 4 i and Ti 2 = t 34 = t 23 /7. It was executed 
for frequency ranging from 0.025 to 1 Hz by varying t 23 between 
17.5 and 0.438 s. The high electric field E H varied from 0.75 to 
0.95 MV/m while T H varied from 60 to 217 °C, and the uniaxial 
stress a H applied during processes 2-3 and 3-4 ranged from 0 to 
33.56 MPa. The low electric field E L was fixed at 0.2 MV/m. The cold 
source temperature T c was passively maintained near room tem¬ 
perature around 22 °C and never exceeded 30 °C. In addition, ther¬ 
momechanical cycling was performed without electric field cycling 
for the above conditions on one of the samples in order to calibrate 
the sample temperature oscillations for different frequencies and 
hot source temperature T H . To do so, a type-K thermocouple was 
bonded with OMEGABOND* 101 to the center of one of the 
3x5 mm 2 faces. 

4.3.4. Olsen cycle 

The Olsen cycle was performed to achieve the maximum energy 
density with cold source temperature T c set at 22 °C and hot 
source temperature T H ranging from 80 to 170 °C. The electric 
fields E l and E H were fixed at 0.2 and 0.75 MV/m, respectively. 
The cycle frequency / was dependent on T H based on the amount 
of time necessary for the sample’s electric displacement to reach 
a minimum or maximum during the heating and cooling processes, 
respectively. Note that the maximum energy density for [001] 
PMN-28PT undergoing the Olsen cycle has previously been 
achieved with sample temperature T cM = 90 °C [42], For both 
cycles, the energy density generated per cycle N D was evaluated 
by numerically integrating experimental data for D vs E according 
to Eq. (2) using the trapezoidal rule. 


“dipping experiments” and the experimental setup described in 

Refs. [31,32], 5- Results and discussion 


4.3.1. Isothermal D-E loops 

Isothermal bipolar D-E hysteresis loops were collected on the 
samples for temperature ranging between 22 and 170 °C and com¬ 
pressive stress varying from 0 to 25.13 MPa. These measurements 
were taken by applying a triangular voltage signal at 0.1 Hz across 
the single crystal samples. The amplitude of the applied voltage 
corresponded to an electric field varying from -0.75 to 0.75 MV/m. 
All measurements were repeated three times to assess repeatability 
and to estimate the experimental uncertainty. 


4.3.2. Differential scanning calorimetry 

Differential scanning calorimetry (DSC) measurements were 
performed using a Diamond DSC (by Perkin Elmer, USA). The spe¬ 
cific heat c p (T) of the PMN-28PT sample was estimated as 


c P (T) 


Qdsc 

Tnf 


( 6 ) 


where Qdsc is the measured heat transfer rate (in W) to achieve con¬ 
stant heating or cooling rate of t = dT/dt (in °C/s) for a sample of 
mass m (in kg). 

The DSC instrument was successfully calibrated using an 
indium standard. The melting temperature at atmospheric pres¬ 
sure and specific phase change enthalpy were measured to be 
160.8 °C and 26.9 J/g, respectively. These values fall within 2.7% 
and 5.6% of the properties reported in the literature [49]. The pro¬ 
cedure to measure the specific heat c p was also validated using an 


5.1. Isothermal bipolar D-E loops 

Fig. 3 plots the isothermal bipolar D-E loops measured at 0.1 Hz 
at temperature (a) 22 °C, (b) 80 °C, (c) 140 °C, and (d) 150 °C and 
under mechanical loading varying between 0 and 25.13 MPa. It 
illustrates the effect of compressive stress on the D-E loops. The 
isothermal D-E loops were closed and consecutive D-E loops over¬ 
lapped for any temperature and compressive stress considered. 
This indicates that leakage current through the sample was negli¬ 
gibly small. The non-linearity in D-E loops, observed as the electric 
field was reduced from 0.75 to 0.0 MV/m under zero stress, corre¬ 
sponded to electric field induced phase transitions [28 . According 
to the E-T phase diagram at 10 Hz [22], tetragonal to monoclinic 
M c phase transition occurs at 80 °C and 0.4 MV/m. Similarly, the 
tetragonal to cubic phase transition at temperature 140, 150, and 
160 °C occurred at electric fields 0.1, 0.18, and 0.25 MV/m, respec¬ 
tively [22]. At T = 170 °C, the sample remained in the paraelectric 
pseudocubic phase for all electric fields considered and compres¬ 
sive stress had negligible effect on the D-E loops (not shown). 

5.2. Specific heat 

Fig. 4 plots the specific heat c p of [001 [-poled PMN-28PT as a 
function of temperature T during heating between 15 and 175 °C, 
followed by cooling from 175 to 15 °C. The sample mass and heat¬ 
ing rate were m = 148.3 mg and t = 10 °C/min, respectively. Similar 
results were obtained at heating rate of 5 °C/min. The measured 
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Electric field, E (MV/m) 


Fig. 3. Isothermal bipolar D-E loops of PMN-28PT Sample 1 at 0.1 Hz for temperatures (a) 22, (b) 80, (c) 140, and (d) 150 °C and compressive stress a between 0 and 
25.13 MPa. 


specific heat around room temperature was in good agreement 
with data reported in the literature [24], Tang et al. [24] reported 
the specific heat at room temperature of PMN-xPT single crystal 
to be 2.5 MJ/m 3 K over a wide range of compositions. This corre¬ 
sponded to 319 J/kg I< assuming the density of PMN-28PT to be 
7847 kg/m 3 . In addition, during heating, the PMN-28PT sample 
exhibited peaks in c p at the phase transition temperatures of 85 
and 148 °C corresponding respectively to the monoclinic M A to 
tetragonal (T) and tetragonal to cubic (C) phase transitions as sug¬ 
gested by the field cooling phase diagram for [001 [-poled PMN- 
28PT [21 j. During cooling, the PMN-28PT sample underwent the 
reverse phase transition sequence with the phase transitions 
occurring at 132 and 69 °C. Such thermal hysteresis [10] has also 
been previously observed in PMN-xPT compositions with x ranging 
from 0% to 25% [25], 


5.3. Sample temperature calibration 

Fig. 5a-d show the sample temperature oscillations as a func¬ 
tion of time for different heater temperatures T H varying from 80 
to 200 °C measured at 0.025, 0.125, 0.5, and 1 Hz, respectively. 
For all heater temperatures, the sample temperature took around 
40 s to reach oscillatory steady-state. For a given frequency, the 
peak to peak temperature span and the maximum and minimum 
sample temperatures T hot and T coM increased with increasing heater 
temperature T H . Fig. 6a and b show the sample temperature cali¬ 
bration for sample temperatures T/, m and T coM as a function of hea¬ 
ter temperature T H for frequency varying from 0.025 to 1 Hz. Each 



Sample temperature, T (°C) 


Fig. 4. Measured specific heat c p of PMN-28PT as a function of temperature 
between 15 and 175 °C. 


data point corresponds to the arithmetic mean of T hot and T coU 
measured over five consecutive cycles in the oscillatory steady- 
state regime. The solid lines correspond to the linear fit of T hot or 
Tcoid vs T h (in °C) for a given frequency. This fit was used to com¬ 
pute T hol or Tcoid for a given heater temperature and cycle fre¬ 
quency. In addition, the sample cold temperature T cou increased 






















84 


I.M. McKinley et al./Applied Energy 126 (2014) 78-89 






Fig. 5. (a) Sample temperature as a function of time for the new cycle performed at (a) 0.025 Hz, (b) 0.125 Hz, (c) 0.5 Hz, and (d) 1 Hz. 


with increasing T H and frequency. In fact, the difference between 
T hot and T coM decreased as the frequency increased to nearly vanish 
at 1 Hz. Furthermore, the sample temperature was found to be 
independent of compressive stress for any given heater tempera¬ 
ture. Note that in the D-E loop experimental setup where the sam¬ 
ple was heated from both top and bottom, its temperature was 
approximately 27 °C less than the heater temperature under 
steady-state conditions due to the low thermal conductivity 
(~0.12 W/m I<) of the electrically insulating Kapton film. 

5.4. Power cycle energy density comparison 

Fig. 7a and b depict, in the D-E diagram, the new thermome¬ 
chanical power cycle performed at 0.025 Hz between cold and 
hot source temperatures T c = 22 °C and T H = 157 °C or T H = 187 °C 
corresponding to sample temperatures (a) T coU = 52 °C and 
T hot = 109 °C or (b) T coM = 59 °C and T h0 , = 129 °C, respectively. For 
both cycles, the electric field was cycled between E L = 0.2 MV/m 
and E h = 0.75 MV/m. A compressive stress of a H = 18.81 MPa was 
applied during processes 2-3 and 3-4. Based on sample tempera¬ 
ture calibration curves (Fig. 6), the sample did not reach the cold 
source temperature of 22 °C, and instead, cooled to approximately 
52 and 59 °C during process 4-1 before process 1-2 was performed. 
Fig. 7 also shows the isothermal bipolar D-E loops previously col¬ 
lected at temperatures near T hot and T coM and under compressive 
stress 0 and <r H . The power cycles shown were vertically translated 


to match the electric displacement of the D-E loop for T hol and a H 
at E h (state 3) as performed by Olsen and Evans [51], The new 
power cycles shown in Fig. 7a and b generated energy densities 
of 24.0 and 33.6 J/L, respectively. The Olsen cycle performed under 
similar operating temperatures and electric fields would yield 
much smaller energy densities considering the area bounded by 
the D-E loops at temperatures T coM and T hot [42], In fact, it is 
evident from Fig. 7a that for T coW = 52 °C and T hot = 109 °C the Olsen 
cycle would generate no energy. For T coM = 59 °C and T hot = 129 °C, 
the Olsen cycle would generate approximately half of the energy 
density of the new cycle. 

Fig. 8 plots the energy density experimentally generated (i) by 
the new thermomechanical power cycle at frequency 0.025 Hz 
and (ii) by the Olsen cycle for T coW ~90 °C at frequency ranging 
from 0.013 to 0.021 Hz as a function of temperature T h0 , varying 
between 80 and 170 °C. Note that data for the Olsen cycle corre¬ 
sponded to the maximum energy density [42], Each data point rep¬ 
resents an average over four cycles and the error bars correspond 
to one standard deviation or 63% confidence interval. The Olsen 
cycle was performed using the dipping method [31] on the same 
PMN-28PT sample used to perform the new cycle. The error bars 
associated with the Olsen cycle were larger because it was not 
automated, unlike the new cycle. Fig. 8 shows that, for T hot 
between 80 and 150 °C, the energy density generated by the new 
power cycle was larger than that generated by the Olsen cycle. In 
fact, the Olsen cycle did not generate positive energy density for 
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Fig. 6. Calibration curves for maximum and minimum sample temperatures (a) T hot 
and (b) T c „m as functions of heater temperature T H for frequency varying from 0.025 
to 1 Hz. The lines correspond to linear fits of the experimental data for a given /. 


T ho i below 130 °C. However, it generated the largest energy density 
of 85.9 J/L/cycle at T H = T h0[ = 170 °C. Similarly, the energy density 
obtained with the new power cycle increased with increasing 
and reached a maximum of 42.6J/L/cycle for T hot = 150 °C. 


5.5. Effect of compressive stress and temperature 

Fig. 9 shows the generated energy density of the new thermo¬ 
mechanical power cycle as a function of applied compressive stress 
a H for heater temperatures T H varying between 60 and 170 °C. 
Here, the frequency was fixed at 0.125 Hz while the electric field 
was cycled between E L = 0.2 MV/m and E H = 0.75 MV/m. It is evi¬ 
dent that the energy density increased nearly linearly with increas¬ 
ing a H for any given heater temperature T H . it also increased 
slightly with increasing heater temperature up to T H = 160 °C. For 
this heater temperature and frequency, the sample temperature 
T co u was 90 °C, corresponding to the temperature at which [001] 
PMN-28PT was in the tetragonal phase and exhibited the highest 
saturation polarization under zero stress [42], On the other hand, 
for T h = 170 °C and /= 0.125 Hz, T coM was 96 °C, corresponding to 
the pseudocubic phase featuring a smaller saturation polarization 
than the T-phase under zero stress. Therefore, the energy density 




Fig. 7. Isothermal bipolar D-E loops and experimental power cycles for PMN-28PT 
for electric field between E L = 0.2 MV/m and E H = 0.75 MV/m and temperatures 
between (a) T cM = 52 °C and T tat = 109 °C or (b) T c0 , d = 59 °C and T ho , = 129 °C. The 
power cycles were performed at frequency f - 0.025 Hz. State 3 of the power cycles 
were vertically displaced to match the corresponding D-E loop at T hot and E H . The 
compressive stress applied during processes 2-3 and 3-4 was equal to 18.81 MPa. 


generated at 0.125 Hz for a given compressive stress was the larg¬ 
est for T h = 160 °C. 

5.6. Effect of frequency 

Fig. 10a shows the energy density generated using the new 
power cycle as a function of heater temperature T H for frequency 
0.125, 0.5, and 1 Hz. The applied compressive stress <j h was fixed 
at 25.13 MPa, the low and high electric fields E L and E H were 
respectively set at 0.2 and 0.75 MV/m, and the cold source temper¬ 
ature T c was 22 °C. it is evident that the energy density increased 
with increasing heater temperature for any given frequency. In 
addition, it decreased only slightly with increasing cycle frequency. 
This was due to the fact that the change in electric displacement as 
a result of a change in compressive stress (process 2-3) occurred 
quickly and was nearly similar for all frequencies. Furthermore, 
longer heating and cooling durations t 2 3 and r 4 i at lower cycle fre¬ 
quency allowed the slow thermal relaxation processes to take 
place. This enabled more charge to build up at the electrode sur¬ 
faces during process 4-1 resulting in the largest N 0 generated for 
T h above 110 °C at 0.125 Hz. In fact, the difference between T coU 
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Fig. 8. Experimentally measured energy density generated with PMN-28PT as a 
function of hot temperature T hot for maximum energy density Olsen cycles and new 
power cycles at 0.025 Hz. The Olsen cycle was performed with zero compressive 
stress ((7 h = 0 MPa) and the new power cycle was performed with a H = 18.81 MPa. 
For all cycles T c = 22 °C, E L = 0.2 MV/m, and E H = 0.75 MV/m. 



Fig. 9. Experimentally measured energy density generated by performing the new 
thermomechanical power cycle on PMN-28PT as a function of compressive stress. 
The heater temperature T H varied from 60 to 170 °C while T c = 22 °C,/= 0.125 Hz, 
E l = 0.2 MV/m, and E H = 0.75 MV/m. 


and T hot was around 10 °C at 0.125 Hz, while it was less than 1 °C at 
higher cycle frequency. 

Fig. 10b shows the power density generated by the new power 
cycle as a function of heater temperatures T H corresponding to the 
data shown in Fig. 10a. It is evident that increasing the cycle fre¬ 
quency resulted in significantly larger power density. In fact, the 
power density P 0 (= N D f) at 1 Hz was nearly ten times larger than 
that at 0.125 Hz for heater temperature T H above 100 °C. This can 
be attributed to the fact that N 0 did not decrease significantly as 
the frequency increased from 0.125 to 1 Hz (Fig. 10a). 

Fig. 11a and b show typical electric displacement changes D-D 4 
vs time f for the new thermomechanical power cycle performed at 
frequency 0.125 and 1 Hz, respectively. For illustration purposes, 
state 4 of the cycle was used as a reference. Fig. 11a and b also 
show the states 1 to 4 of the cycle corresponding to those shown 
in Fig. lc. In both cases, the heater temperature T H was maintained 




Fig. 10. Experimentally measured (a) energy density and (b) power density 
generated by performing the new power cycle on PMN-28PT as a function of 
heater temperature T H for frequency varying from 0.125 to 1 Hz with T c = 22 °C, 
E l = 0.2 MV/m, E„ = 0.95 MV/m, and <7„ = 25.13 MPa. 


at 160 °C, the low and high electric fields E L and E H were respec¬ 
tively 0.2 and 0.75 MV/m, and the compressive stress was 
tr H = 25.13 MPa. Fig. 11a clearly illustrates the piezoelectric and 
pyroelectric contributions to the cycle at 0.125 Hz. During process 
4-1, approximately 85% of the rise in electric displacement from 0 
to 6 pC/cm 2 occurred in the first 0.25 s. This change can be attrib¬ 
uted to the change in compressive stress from a H to 0 MPa. Then, 
the remaining 15% of the change in electric displacement during 
process 4-1 occurred between 0.25 and 3.5 s and can be attributed 
to cooling of the sample from T hot = 101 °C to T coM = 90 °C. In addi¬ 
tion, Fig. lib shows a similar response to the reduction in com¬ 
pressive stress during process 4-1 as the electric displacement 
increased from 0 to 6 pC/cm 2 in the first 0.25 s. However, at 1 Hz 
the sample did not have time to experience cooling during process 
4-1 and no additional changes in electric displacement took place. 
As a result, the generated energy density as well as the overall 
change in electric displacement between states 4 and 2 for the 
cycle at 1 Hz were approximately 85% of those at 0.125 Hz. This 
indicates that the piezoelectric contribution to the new cycle was 
independent of cycle frequency between 0.125 and 1 Hz while 
the pyroelectric contribution decreased with increasing frequency. 
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Fig. 11. Experimentally measured electric displacement vs time for the new power 
cycle for frequency (a) 0.125 and (b) 1 Hz. The temperature T H was maintained at 
160 °C, the low and high electric fields E L and E H were set at 0.2 and 0.75 MV/m, 
while <t h was 25.13 MPa. 

5.7. Maximum power density 

Fig. 12 shows five consecutive new cycles performed on 
PMN-28PT in the D-E diagram corresponding to the maximum 
generated power density. The heat sink and heat source tempera¬ 
tures T c and T H were set at 22 and 130 °C, respectively. The 
high compressive stress a H was 33.56 MPa. The low and high 
electric fields were 0.2 and 0.95 MV/m, respectively. The cycle 
frequency was 1 Hz and the resulting sample temperature was 
ThotceTcoid = 85 °C, corresponding to the phase boundary between 
the monoclinic phase and the highly polarized tetragonal phase. 
The consecutive cycles nearly overlapped indicating that the cycle 
was highly repeatable resulting in power density of 41.3 ± 0.4 W/L. 

5.8. Material efficiency 

Table 1 summarizes the operating conditions, energy inputs Q ln 
and W in , generated energy density N D , and the material efficiency ;/ 
obtained experimentally for various conditions of the new power 
cycle performed on PMN-28PT. It also gives conditions correspond¬ 
ing to the maximum material efficiency of the Olsen cycle. The 


Fig. 12. D-E path of five consecutive cycles of the highest experimentally measured 
power density generated with PMN-28PT with T c = 22°C, T H = 130°C, /= 1 Hz, 
E l = 0.2 MV/m, E h = 0.95 MV/m, and a„ = 33.56 MPa. 

material efficiency of the thermomechanical power cycle and of 
the Olsen cycle was estimated using Eq. (3). The thermal energy 
Q,„ was estimated according to Eq. (4) using the measured density 
and specific heat c p presented in Fig. 4. The mechanical energy 
input W in was computed from Eq. (5) using the area enclosed by 
the stress-strain curves at 85 °C measured at 0.2 and 0.95 MV/m 
for [001] PMN-28PT single crystals [52], When the new cycle 
was performed at low frequency, the large heat input resulted in 
small efficiencies. These results illustrate the importance of operat¬ 
ing under small temperature swings and heat input to improve the 
cycle efficiency. 

Moreover, for a given temperature swing, the efficiency was 
highly dependent on the operating temperatures. For example, 
the two cycle performed at 0.5 Hz in Table 1 both had similar 
energy density N D , the same W,„, and a temperature swing 
Tha - Tcoid = 1 °C, but their efficiencies differed by nearly a factor 
of two. This difference can be attributed to the operating tempera¬ 
tures. One cycle was performed between 82.3 and 83.2 °C where 
c p (T) exhibited a phase transition peak and a large thermal hyster¬ 
esis (Fig. 4) while the other was performed between 99.9 and 
100.9 °C, where no phase transition occurred. This emphasizes 
the importance of choosing the operating temperatures such that 
the heat input is minimized and the energy density is maximized 
in order to maximize the material efficiency >/. This can be achieved 
by operating under adiabatic conditions at a bias temperature near 
the Curie temperature [52], In fact, the maximum power density 
cycle performed at 1 Hz (Fig. 12) achieved an efficiency larger than 
64%. In this case, the heat input was Q m = 9 kj/m 3 and the mechan¬ 
ical work performed was W in = 55 kj/m 3 . Although the sample 
underwent small temperature fluctuations T hot - T roW ~0.1 °C dur¬ 
ing the cycle, the corresponding thermal energy input Q in was 
comparable to the mechanical energy W in . 

In addition, at low frequencies, it is important to choose T cou 
such that heating and applying compressive stress have a similar 
effect on the electric displacement. The two cycles performed at 
0.125 Hz in Table 1 illustrate this fact. The cycle performed 
between 90.4 and 103.1 °C had a material efficiency of 9.4%. In this 
case, the material was in the highly polarized tetragonal phase at 
Tcoid = 90.4 °C and E H = 0.75 MV/m [22], Then, both heating and 
the addition of compressive stress during process 2-3 decreased 
the electric displacement of the material. On the other hand, for 
the cycle performed between 74.4 and 85.1 °C, the sample was 
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Table 1 

List of operating conditions, energy inputs Q, n and W m , generated energy density N D , and material efficiency p for the new cycle and the maximum Olsen cycle efficiency. In all 
cases, the low electric field was E L = 0.2 MV/m. 


Cycle type 

/(Hz) 

E h (MV/m) 

r„(°c) 

That (°C) 

Gom (°C) 

<Jh (MPa) 

w,„ (kj/m 3 ) 

Qtn (kj/m 3 ) 

N d (kj/m 3 ) 

»/(%) 

Olsen 

0.021 

0.75 

- 

170 

90 

0 

0 

540 

86 

15.9 

TM 

0.025 

0.75 

157 

108.5 

51.9 

18.81 

30 

1202 

24 

2.0 

TM 

0.025 

0.75 

107 

73.9 

40.2 

18.81 

30 

292 

15 

4.7 

TM 

0.125 

0.75 

160 

103.1 

90.4 

25.13 

41 

301 

32 

9.4 

TM 

0.125 

0.75 

130 

85.1 

74.4 

25.13 

41 

469 

31 

6.1 

TM 

0.5 

0.75 

160 

100.9 

99.9 

25.13 

41 

21 

27 

43.5 

TM 

0.5 

0.75 

130 

83.2 

82.2 

25.13 

41 

62 

26 

25.2 

TM 

1 

0.75 

160 

104.4 

104.3 

25.13 

41 

2 

26 

60.5 

TM 

1 

0.75 

130 

85.4 

85.3 

25.13 

41 

9 

25 

50.0 

TM 

1 

0.95 

130 

85.4 

85.3 

33.56 

55 

9 

41 

64.1 


near the phase boundary between monoclinic M c and tetragonal at 
T c oid = 74.4 °C and E H = 0.75 MV/m [22], Then, during process 2-3, 
the electric displacement decreased with applied compressive 
stress and increased with heating. In this case, the material effi¬ 
ciency was 6.1%. Both of these cycles yielded a similar energy den¬ 
sity, however one consumed an additional 168kJ/m 3 of thermal 
energy per cycle. 

Finally, note that the efficiency reported corresponds to the 
material’s ability to convert thermomechanical energy into elec¬ 
tricity. The efficiency of a device implementing the new cycle on 
PMN-28PT is expected to be significantly lower. In addition, W,„ 
was estimated based on properties of PMN-28PT measured at 
85 °C [52], In practice, the mechanical work was performed at 
two different temperatures and the stress-strain behavior may 
have deviated from that at 85 °C [52]. Moreover, Q,„ was estimated 
based on the specific heat c p (T) measured under zero mechanical 
stress. In practice, the heat input occurred under large uniaxial 
compressive stress known to affect the phase transition tempera¬ 
tures and the corresponding peaks in c p (T ). Furthermore, note that 
the material efficiency of the Olsen cycle reached 48% of the Carnot 
efficiency between T c = 22 °C and T H = 170 °C. 

6. Conclusion 

This paper presented the concept and experimental implemen¬ 
tation of a novel thermomechanical power cycle performed on 
[001 [-poled PMN-28PT single crystals. The cycle consisted of four 
processes, namely (i) increasing the applied electric field under 
constant temperature T cM and zero applied stress, (ii) simulta¬ 
neously heating up to temperature T hot and compressing the mate¬ 
rial under compressive stress a H under large electric field, (iii) 
decreasing the applied electric field under constant temperature 
Thot and compressive stress a H , and finally (iv) simultaneously 
removing the applied stress and cooling the material down to 
Tcom at small electric field. Maximum energy and power densities 
of 41 J/L/cycle and 41 W/L were achieved with the cycle performed 
at 1 Hz for heat sink and heat source temperatures of 22 °C and 
130 °C and electric field cycled between 0.2 and 0.95 MV/m under 
compressive stress 33.56 MPa. These conditions also yielded the 
maximum material efficiency of converting thermomechanical 
energy into electric energy of 64.1%. The energy density of the 
new cycle was found to be nearly independent of cycle frequency, 
while the power density and efficiency increased with increasing 
cycle frequency. In addition, both the energy and power densities 
increased with increasing compressive stress. Furthermore, the 
cycle was able to convert thermal and mechanical energies into 
electrical energy at all temperatures considered. It is also capable 
of maintaining constant power output for changing hot source 
temperature by adjusting the applied stress. Future work should 
consist of implementing this cycle into autonomous and integrated 
devices. 
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Glossary 

Coercive electric field: is the electric field applied to a pyroelectric material corre¬ 
sponding to zero electric displacement. 

Curie temperature: is the temperature at which a ferroelectric material undergoes a 
phase transition from ferroelectric to paraelectric. 

D-E loops: are electric displacement versus electric field curves of pyroelectric 
materials that indicate material properties such as, large-field permittivity, 
coercive electric field, and remnant and saturation polarizations. 

Electric displacement: is the electric charge per unit electrode area within a dielec¬ 
tric material. It consists of the permanent and the induced electric dipole 
moments in the material. 

Ferroelectric materials: is a subclass of pyroelectric materials having the ability to 
switch the direction and magnitude of the spontaneous polarization by apply¬ 
ing an electric field above the coercive electric field. 

Morphotropic phase boundary: refers to phase transitions due to changes in com¬ 
position of ferroelectric materials. It usually corresponds to the phase transition 
between the tetragonal and rhombohedral ferroelectric phases. 

Olsen cycle: is a cycle that is performed on a pyroelectric material to utilize time- 
dependent temperature oscillations to convert thermal energy directly into 
electricity. It is sometimes known as the Ericsson cycle. 

Large-field permittivity: is the slope of the D-E loop of a pyroelectric material at 
large electric field divided by the vacuum permittivity e 0 = 8.854 x 10 -12 F/m. 

Leakage current: refers to the transport of charges accumulated at the surface of a 
pyroelectric material through its body. 

Paraelectric: refers to a state of a pyroelectric material that is in a crystal phase in 
which the electric dipoles are unaligned resulting in zero spontaneous 
polarization. 

Piezoelectric materials: are dielectric materials in which the electric charge in the 
material and the mechanical deformation are related, i.e., an applied electric 
field induces a mechanical deformation and a mechanical deformation results 
in a change in electric displacement. 

Pyroelectric materials: is a subclass of piezoelectric materials that possess a tem¬ 
perature-dependent spontaneous polarization. 

Remnant polarization: is the polarization exhibited by a pyroelectric material under 
zero applied electric field. 

Saturation polarization: is equal to the electric displacement in the linear fit of the 
D-E loops of a pyroelectric material at large field extrapolated to zero electric 
field. 

Spontaneous polarization: is the average electric dipole moment per unit volume of 
a pyroelectric material in the absence of an electric field. 

Thermal hysteresis: refers to the difference in phase transition temperatures of a 
pyroelectric material during heating and cooling. 



